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ABSTRACT 

We have observed the recently discovered rich star cluster GLIMPSE-C01 for 46 ks with the Chandra X-ray 
Observatory. Seventeen X-ray sources with luminosities > 0.6 x 10 31 ergs s _I were discovered, one of which 
is likely a quiescent low-mass X-ray binary. The spatial distribution of these X-ray sources with respect to 
the NIR and IR images of the cluster, combined with the luminosity and spectral hardness distributions of the 
sources, provide strong evidence that GLIMPSE-CO 1 is a rich Galactic globular cluster. 
Subject headings: globular clusters: individual (GLIMPSE-C01) — X-rays: binaries 



1. INTRODUCTION 

The star cluster "GLIMPSE-C01" was discovered only re- 
cently with the Spitzer Space Observatory during the Galactic 
Legacy Infrared Mid Plane Survey (Kobulnicky et al. 2005), 
and by cross-correlating ASCA Galactic-ridge X-ray sources 
with 2MASS images (Simpson & Cotera 2004). This heavily 
obscured object lies directly in the Galactic plane and is very 
likely a new rich Galactic globular cluster (see Fig.|TJ>. The 
cluster center has Galactic coordinates / = 31° and b = —0.1°, 
placing it within 10 pc of the Galactic midplane, as close as 
any cluster listed in the Harris (1996) catalog 4 . The visual ex- 
tinction is estimated to be Ay — 12- 18 magnitudes (Kobul- 
nicky et al. 2005), which corresponds to an absorbing column 
density of 2.7 ±0.5 x 10 22 cm" 2 (Predehl & Schmitt 1995). 

Kobulnicky et al. (2005) found that GLIMPSE-C01 is rela- 
tively close to the Earth at a distance of 3-5 kpc, has a mass of 
~ 1-3 x 1O 5 M and possibly higher, and seems to be highly 
centrally condensed, with a half-light radius of ^0.5-1 pc and 
a core radius that is a factor of ~4 times smaller (see also 
Ivanov, Kurtev, & Borissova 2005). These half-light and core 
radii are smaller than those for many of the interesting clusters 
that have previously been studied with Chandra, implying a 
very high density in the cluster core and a concomitantly high 
production rate for X-ray sources. 

The possibility that GLIMPSE-CO 1 might be either an old 
open cluster or a more distant "super star cluster" comprised 
of many young, massive stars is extensively discussed by 
Kobulnicky et al. (2005). However, the authors rather con- 
vincingly discount these alternatives on the basis that (i) no 
radio emission or NIR emission lines are observed, (ii) the 
cluster has a well-populated giant branch, and (iii) it is lo- 
cated w 6 kpc from the Galactic center and old open clus- 
ters are extremely scarce at galactocentric radii of < 7.5 kpc 
(Friel 1995), probably because of the disruptive effects of gi- 
ant molecular clouds. GLIMPSE-C01 is the richest cluster 
uncovered thus far in the Spitzer midplane survey covering 92 
square degrees. 

The field near GLIMPSE-C01 has not previously been ob- 
served by Chandra but was observed with XMM-Newton for 
15 ksec (P. Woods, private communication). In addition, the 
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ASCA survey of the Galactic Ridge (Sugizaki et al. 2001) lists 
source AX J 184848-0219 within 45" of the GLIMPSE-C01 
cluster center. Neither the XMM nor ASCA observation had 
the angular resolution required to resolve the X-ray sources in 
the core of this cluster. The XMM image reveals extended 
emission at the location of the globular cluster that appar- 
ently represents the collection of point sources identified in 
this work. 

Here we report the results of a 46 ks Chandra observation 
of the GLIMPSE-CO 1 cluster. Seventeen X-ray sources with 
luminosities above ^0.6 x 10 31 ergs s" 1 were discovered. We 
show that their spatial and spectral properties indicate that 
GLIMPSE-CO 1 is a globular cluster as opposed to an open 
cluster. 

2. CHANDRA OBSERVATIONS 

GLIMPSE-CO 1 was observed for 46 ks on 2006 August 15- 
16 (ObsID 6587) with the Advanced CCD Imaging Spectrom- 
eter (ACIS) on Chandra with the telescope aim point on the 
back-side illuminated S3 chip. The data were taken in timed- 
exposure mode with the standard integration time of 3.24 s 
per frame and telemetered to the ground in faint mode. 

Data reduction was performed using the CIAO 3.4.1 soft- 
ware provided by the Chandra X-ray Center 5 . The data were 
reprocessed using the CALDB 3.4.0 set of calibration files 
(gain maps, quantum efficiency, quantum efficiency unifor- 
mity, effective area) including a new bad-pixel list made with 
the acis_mn_hotpix tool. The reprocessing was done without 
including the pixel randomization that is added during stan- 
dard processing. This omission slightly improves the point 
spread function. The data were filtered using the standard 
ASCA grades (0, 2, 3, 4, and 6) and excluding both bad pixels 
and software-flagged cosmic ray events. Intervals of strong 
background flaring were searched for, but none were found. 
The extraction of counts and spectra and the generation of re- 
sponse files and background spectra were accomplished with 
ACIS Extract v3.131 (Broos et al. 2002), which calls many 
standard CIAO routines. All extraction regions were con- 
structed to match the shape of the point-spread function (PSF) 
and enclose 90% of the PSF. 

3. X-RAY IMAGE AND SOURCE PROPERTIES 

The CIAO wavelet-based wavdetect tool was employed for 
source detection in the 0.5-6.0, 0.3-10.0, and 2-8 keV bands. 
We detected 42 sources on the entire S3 chip. Of these 
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FIG. 1 . — 2' X 2' images of the field of GLIMPSE-C01 with the 36" radius half-light region indicated by the white circle. Left panel: Adaptively smoothed 
Chandra image with sources indicated. Right panel: True-color image made from the J, H, and K bands of an NTT image (adapted from Ivanov et al. 2005). 



sources, 12 lie within the 36" half-light radius of GLIMPSE- 
C01 . From the density of sources outside the half-light radius, 
we expect 0-1 non-cluster sources inside the half-light radius. 
We also examined adaptively-smoothed images (made from 
the CIAO csmooth tool) for possible point sources. We found 
five additional candidate sources this way, but some of these 
may be spurious. The full source list is given in Table 1, with 
sources numbered according to detected counts in the 0.5-8 
keV range. 

We fit all sources with absorbed power-law spectral models 
in Sherpa (Freeman et al. 2001) using Cash (1979) statistics, 
which are appropriate for low count data. We fixed the col- 
umn density to 2.7 x 10 22 cm" 2 . For each source, we simul- 
taneously fit the source spectrum and associated background 
spectrum, which was constructed from a source-free annulus 
around the source. We performed the simultaneous fitting 
because Cash statistics are not appropriate for background- 
subtracted spectra. The best fit power-law photon index and 
unabsorbed 0.5-8 keV luminosity, assuming a distance of 3.7 
kpc (Ivanov et al. 2005) are given in Table 1. 

The X-ray image of GLIMPSE-C01 is shown in Fig[T] 
along side an NTT image constructed from the J, H, and K 
bands (courtesy of V. Ivanov & J. Borissova, private commu- 
nication; adapted from Ivanov et al. 2005) Both images are 
color coded for photon energy and intensity. The X-ray source 
designations are shown in the Chandra image. 

Source positions, luminosities, and spectral properties are 
summarized in Table 1. In addition, we list the results of 
two statistical tests performed by ACIS Extract. The first 
is the Kolmogorov-Smirnov probability that the source spec- 
trum and background spectrum are drawn from the same par- 
ent distribution; values close to unity indicate they are. We 
note that because a background spectrum is constructed from 
a local annulus around a source and because there is diffuse 
emission which may be unresolved point sources, a similarity 
of the source and background spectral shapes is not necessar- 
ily sufficient grounds to regard a particular source as spurious. 



The second statistic in the table is the Kolmogorov-Smirnov 
probability that the source's photon arrival times are consis- 
tent with a constant-intensity source; values close to zero in- 
dicate a variable source. 

The sources are centrally concentrated around the clus- 
ter core, the radius of which we estimate to be ~9" (see 
§4.2). This is quite similar to the spatial distribution of 
low-luminosity X-ray sources found in other Galactic glob- 
ular clusters (e.g., Pooley et al. 2003; Heinke et al. 2005). 
Overall, the range of spectral luminosities, hardness ratios, 
and spatial distribution of the sources clearly associated with 
GLIMPSE-CO 1 are highly reminiscent of a Galactic globular 
cluster rather than those of an open cluster (see, e.g., Fig. 1 1 
of Kobulnicky et al. 2005; van den Berg et al. 2004). 

In addition to the point sources there is also a significant 
component of diffuse emission nearly aligned with the clus- 
ter center and very likely directly associated with the cluster. 
After subtracting out the contribution from the detected point 
sources and the background, we find an excess of 513 ±45 
counts [0.5-8 keV] within the optical half-light radius. This 
is comparable to (only about 15% greater than) the total num- 
ber of counts from all detected point sources (the sum of the 
counts listed in Table 1, corrected for the 90% PSF extrac- 
tion regions). In the 0.5-2 keV band, the diffuse emission is 
82 ±20 counts, and, in the 2-8 keV band, it is 432 ±40 counts. 
The diffuse X-ray emission has an asymmetric morphology 
that appears to peak near the locus of X7, X5, X6, and XI 1, 
which is offset about 8" from the optical center of the clus- 
ter. This may be due to differential extinction across the field 
(Kobulnicky et al. 2005; Ivanov et al. 2005); however, even 
in the 3-8 keV band, which is insensitive to extinction, there 
appears to be asymmetry. The X-ray excess has a "half-light" 
radius of about 20". 

4. DISCUSSION 
4.1. X-ray Sources 
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TABLE 1 
X-ray Sources in GLIMPSE-C01 



Source 


RA (J2000) 


Decl. (J2000) 


Error (") a 


Counts 




PL Index 


KSbkgd d 


IV » 3 eonstaoi 


XI 


18:48:50.261 


-01:29:28.57 


0.07 


68 


9.86 


0.1 


5.7x 10 


0.70 


X2 


18:48:48.867 


-01:29:59.56 


0.07 


56 


63.0 


4.6 


1.3x 10~ 12 


0.85 


X3 


18:48:49.231 


-01:30:00.05 


0.08 


43 


6.21 


0.1 


3.9 x 10~ 3 


0.97 


X4 


18:48:49.799 


-01:29:29.96 


0.09 


37 


6.17 


-0.4 


3.5 x 10~ 3 


0.64 


X5 


18:48:49.949 


-01:29:48.74 


0.09 


33 


4.60 


0.2 


1.7 x 10~' 


0.14 


X6 


1 O AO Af~\ /L O A 

18:48:49.684 


-01:29:56.49 


0.09 


33 


4.26 


0.7 


2.6x 10 -2 


0.71 


X7 


18:48:50.163 


-01:29:53.15 


0.10 


30 


4.29 


1.9 


9.2x 10" 4 


0.06 


X8 


18:48:49.505 


-01:29:40.09 


0.12 


22 


3.20 


-0.1 


6.9 xl0~ 3 


0.95 


X9 


18:48:52.012 


-01:29:39.48 


0.16 


14 


1.73 


1.2 


4.9 x 10" 2 


0.60 


X10 


18:48:49.559 


-01:29:47.51 


0.18 


10 


1.22 


0.5 


2.6x nr 1 


2xl0- 6 


Xll 


18:48:50.126 


-01:30:00.70 


0.18 


10 


1.21 


1.5 


8.7 xl0~ 2 


0.47 


X12 


18:48:49.515 


-01:29:59.14 


0.18 


9 


1.05 


1.0 


2.3x 10-' 


0.29 


X13 


18:48:51.147 


-01:29:44.63 


0.22 


7 


1.04 


-0.1 


3.7x 10"' 


0.19 


X14 


18:48:47.932 


-01:29:43.62 


0.21 


7 


0.51 


10 


9.6x10-' 


0.32 


X15 


18:48:49.243 


-01:30:10.21 


0.22 


6 


0.85 


1.9 


7.2x10-' 


0.69 


X16 


18:48:50.505 


-01:29:36.45 


0.28 


5 


0.53 


1.5 


9.5 x 10"' 


0.05 


X17 


18:48:48.183 


-01:29:58.35 


0.25 


5 


0.57 


1.1 


8.5x10-' 


0.63 



a These statistical uncertainties do not take into account the ~0" 4 systematic uncertainty in Chandra's pointing accuracy. 13 Unabsorbed 0.5-8 keV luminosity in 
units of 10 3 ' ergs s-'. c Best fit power law photon index. d Kolmogorov-Smirnov probability that the source spectrum and background spectrum are drawn from 
the same parent distribution; values close to unity indicate they are. c Kolmogorov-Smirnov probability that the source's photon arrival times are consistent with 
a constant-intensity sources; values close to zero indicate a variable source. 



The 46-ksec Chandra observation has allowed us to take 
a census of the X-ray emitting objects in this cluster (with 
L x > 0.6 x 10 31 ergs s" 1 ). We find 17 X-ray sources within the 
half-light radius of the cluster. From the astrometric compar- 
ison that we have made between the X-ray source positions, 
and the crowded IR and NIR field, we were unable to conclu- 
sively associate any of the bright stellar members of the clus- 
ter with the X-ray sources; there is one preliminary match, but 
given the crowded field this might be a chance superposition. 
Since the bright IR sources are likely stellar giants, this lack 
of a correlation is not unexpected. Based on identifications 
of X-ray sources in this luminosity regime in other globular 
clusters, the counterparts are usually quiescent low-mass X- 
ray binaries (qLMXBs; Hertz & Grindlay 1983, Grindlay et 
al. 2001, Pooley et al. 2003), cataclysmic variables (CVs; Ed- 
monds et al. 2003, Knigge et al. 2002, Pooley et al. 2002), 
millisecond radio pulsars (MSPs; Freire et al. 2003, Ransom 
et al. 2005), or chromospherically active main-sequence bina- 
ries (e.g., BY Draconis systems). 

Above a few times 10 31 erg s" 1 , we would expect only 
qLMXBs and CVs, with the qLMXBs spectrally softer than 
the CVs. Given the spectral softness of X2 and its intrinsic L x , 
we conclude that it is likely a qLMXB. The other sources are 
likely a mixture of CVs, MSPs, and BY Dra systems. Among 
these sources, X10 is notable for its extreme variability. Seven 
of its 10 detected counts were recorded in a brief 200 s inter- 
val, with the other three counts coming in the next 8 ks. The 
luminosity given in Table 1 is the average luminosity through- 
out the 46 ks observation; during the 200 s flare, X10 achieved 
L x ~2 x 10 33 erg s _1 . 

4.2. GLIMPSE-C01 Dynamics 

The 100-fold overabundance, per unit mass, of interact- 
ing binaries in globular clusters has been known for 30 years 
(Clark 1975; Katz 1975), and is related to dynamical encoun- 
ters between collapsed stars (neutron stars, white dwarfs, or 
black holes) and cluster field stars or binaries (e.g., Hut et al. 
1992; Pooley et al. 2003, and references therein). The specific 
types of encounters that have been studied include two-body 



tidal capture (Fabian, Pringle, & Rees 1975; Press & Teukol- 
sky 1977; DiStefano & Rappaport 1992, 1994), three-body 
exchange encounters (e.g., Sigurdsson & Phinney 1995; Ra- 
sio, Pfahl, & Rappaport 2001), and direct collisions (Verbunt 
1987; Ivanova et al. 2005). The total rate, T, of close encoun- 
ters between two types of objects is given by: 

no.1no.2V p\rl 3/22 n , 

I oc oc oc pJ r c , (1) 

g a 

where «o,i an d n o.2 are the central densities for the two types 
of objects, V is the effective volume of the interaction region, 
and a is the stellar velocity dispersion. The second propor- 
tionality assumes that both «o 1 an d n o.2 are proportional to 
the central mass density, po, and that the relevant interaction 
volume is that of the cluster core of radius r c ; the third pro- 
portionality uses the virial relation a 1 oc po^.- 

Thus far, Chandra has been used to observe ^25 Galactic 
globular clusters, with typical exposure times of ~ 20-80 
ks; the cluster 47 Tuc has been observed more extensively. 
X-ray luminosities at the limits of sensitivity are typically a 
few x 10 3() ergs s" 1 . A plot of the observed numbers of X-ray 
sources per unit mass in a given cluster versus the estimated 
value of r per unit mass for that cluster is shown in Fig. [2] 
This shows a strong correlation between the number of X-ray 
sources in a globular cluster and T. 

If we assume that GLIMPSE-CO 1 has a mass in the range 
of M c ] uster ~ (1— 3) X 10 5 Mq we can estimate po as follows. 
From Fig. 10 of Kobulnicky et al. (2005), we estimate that 
the surface brightness falls to ~l/2 its central value at a radial 
distance of 9", and we take this as a good approximation to the 
core radius, r c . We have also constructed a radial light profile 
from the Ivanov et al. (2005) NTT image, fit the region out 
to 8 = 1' with a simple Hubble profile (eq. 2), and find a very 
similar value for r c of ^7" ± 2". We note that this value for r c 
that we have adopted differs significantly from the values for 
this parameter cited by Kobulnicky et al. (2005) and Ivanov et 
al. (2005). We further take the distance to GLIMPSE-C01 to 
be 4± 1 kpc. This leads to an estimate of the physical size of 
the core radius of r c ~ 0. 17 ± 0.04 pc. For the half-light radius 



4 



Pooley et al. 



100 





t 


r 


♦ ] ] 







where we have assumed a V-band mass to light ratio of 2, and 
expressed the central density in terms of L pc~ 3 . 

We now use these values for po and r c , and adopt a plausible 
value of a = 10 km s" 1 , to compute a value of the formation 
rate parameter T of 
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FIG. 2. — Number of globular cluster X-ray sources with L x > 4 X 
10 31 ergs/sec per unit mass, n x , versus the encounter rate per unit mass 7 
of the cluster (adapted from Pooley & Hut 2006). Two dozen clusters are 
represented in this figure, but several of the low-7 clusters are binned up 
to improve statistics. The blue line indicates the relationship for the CV- 
dominated population II of Pooley & Hut. The red point is a preliminary 
location on this diagram for the GLIMPSE-CO 1 cluster, with 7 sources above 
the L x limit, T = 400+5n!!, and M = 3 X 10 5 M . 



(the radius encircling half the light in the projected image), we 
adopt the value cited in Table 1 of Kobulnicky et al., viz, Of, — 
36", or r/, = 0.70 ±0.17 pc. Given the roughness of some of 
these estimates, we feel justified in using the modified Hubble 
profile (e.g., Binney & Tremaine 1987) for the cluster surface 
brightness to relate po with r/, and M c i uster : 



£(/?) = 



2p r c 



[l + (R/r c ) 2 ] 



From this, we find that 



M. 



cluster 



47rr 3 ln[l+(r /l /r c ) 2 ] ' 
If we plug in the values we adopted above, we find 
log po = 5.6 ±0.3 [units of L© pc~ 3 ] , 



(2) 



(3) 



(4) 



with the same normalization used in Pooley et al. (2003), 
namely that T / 1 00 is roughly the number of LMXBs expected 
in the cluster. This value of V, after normalizing to the cluster 
mass, is plotted in Fig. 2, along with its considerable uncer- 
tainties. As can be seen, GLIMPSE-C01 is consistent with the 
general trend of other Galactic globular clusters. 

Given that we have shown that GLIMPSE-CO 1 has a high 
production rate for X-ray sources, it appears that this cluster 
is a good candidate to search for msec radio pulsars (Freire et 
al. 2003; Ransom et al. 2005). 
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